AD-A092  930  PRINCETON  UNIV  N  J  DEPT  OF  CHEMISTRY 

STUDIES  IN  MOLECULAR  COLLISION  PHENOMENA, <U> 
I960  H  RABITZ 
UNCLASSIFIED  PR-9 


NOOO14-75-C-O470 

NL 


DDC  ffl-E  COPtj  AD  A0929  3  0 


lEVO? 


Princeton  University 
Department  of  Chemistry 


n  s 


Progress  Report 


under  Contract  N00014-75-C-0478 


Sponsored  by 


Office  of  Naval  Research 


STUDIES  IN  MOLECULAR  COLLISION  PHENOMENA 


for  the  period 

November  1,  1979  -  December  5,  1980 


Herschel  Rabitz 
Principal  Investigator 


ll 


DT1C 

ELECTE 
DEC!  awo 


80  1*  12  0*1 


_ -  s*  v\ 


-  -■ 


SECURITY  CLASSIFICATION  OF  THIS  PACE  'MTien  Data  Fntarad) 


m 

/ 


i  *  ' 


REPORT  DOCUMENTATION  PAGE 


L  RE-PORT  NUMBER 


PR-9 


JL. 


2.  GOVT  ACCESSION  NO 


vUr/im 


4^  TITLE  (and  Subtitle) 

Studies  in  Molecular  Collision  Phenomena  $ 


7.  authors; 

Herschey  Rabitz 


0 "  )Jt  7  fcj  NOA&Vh  V  ^ 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 


Princeton  University 
Princeton,  New  Jersey  08544 


It.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 


+.  1  /V.w  17-5  D* 


14.  MONITORING  AGENCY  NAME  &  ADDRESS^//  different  from  Controlling  Office) 


Office  of  Naval  Research 
Arlington,  Virginia  22217 
Attention:  Code  421 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.  RECIPIENT'S  CAT  ALOG  NUMBER 


Ha 


S.*  TYPE  OF  REPORT  6  PERIOD  COVEREO 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GnANT  NUMBER^; 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  A  WORK  UNIT  NUMBERS 


II)  HJ-. 


/ 


12.  REPORT  DATE 


13.  NUMBER  OF  PAGES 


15.  SECURITY  CLASS,  (of  thie  report) 


15a.  DECLASSIFICATION  -  DOWNGRADING 
SCHEOUt-E 


G.  OSSTRIOUTlOH  ST  AT  EM  t  N  T  < of  this  Report ) 

Approved  for  public  release:  distribution  unlimited 


17.  DISTRIBUTION  STATEMENT  (of  the  ebatrect  entered  In  Block  20,  if  different  from  Report) 


ie.  supplementary  notes 


19.  KEY  WORDS  (Continue  on  reverse  aide  It  neceeeary  and  Identify  by  block  number) 

molecular  collision  theory,  relaxation  processes,  energy  transfer, 
chemical  kinetics 


20.  ABSTRACT  (Continue  on  reverse  side  If  neceeeary  and  Identify  by  block  number) 

Molecular  dynamics  studies  at  Princeton  are  reviewed.  The  areas  of 
research  involved  are.' 

1.  molecular  collisions 

2.  radiation-molecule  interaction 

3.  sensitivity  theory 

4.  finite  elements  in  quantum  mechanics  • 


DD  |  JAN *73  1473  EDITION  OF  I  NOV  85  IS  OBSOLETE 
S/N  010 2* 014s  8601  | 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  Data  Bntatad) 


j  .j  jt  ^ 


J.P 


1 


I.  Contract  Description 

The  research  is  concerned  with  molecular  interactions  and  how  they 
manifest  themselves  in  a  variety  of  physical  phenomena,,  Emphasis  is  on 
attaining  a  better  fundamental  understanding  of  molecular  dynamics  as  well 
as  on  the  study  of  practical  applications. 

II.  Scientific  Problem 

The  aim  of  this  research  is  to  achieve  a  detailed  understanding  of 
the  mechanism  of  chemical  reaction  dynamics  and  energy  transfer  at  a  mole¬ 
cular  level.  This  information  is  of  fundamental  importance  to  molecular 
physics  and  chemistry  as  well  as  being  essential  for  engineering  design 
studies  of  many  practical  devices.  The  research  is  focussed  on  discerning 
the  fate  of  excited  molecules,  particularly  in  a  gaseous  medium.  Such 
excited  molecules  could  react  or  transfer  their  excitation  to  other  mole¬ 
cules  upon  collisional  impact.  In  addition,  an  excited  polyatomic  molecule 
could  undergo  self-relaxation  or  scrambling  of  its  internal  energy  even 
without  a  collision.  The  basic  computed  quantities  are  rate  constants 
for  these  processes.  An  effort  is  under  way  to  map  out  the  relationship 
between  observable  laboratory  behavior  and  the  underlying  input  forces 
between  the  molecules.  In  addition  to  achieving  a  basic  theoretical 
understanding  of  these  processes,  the  research  is  being  used  to  provide 
a  framework  for  the  analysis  and  inversion  of  appropriate  measurements 
to  yield  detailed  molecular  rate  information. 
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Scientific  and  Technical  Approach 


The  varied  nature  of  the  research  has  necessitated  a  multi-faceted 
approach.  For  the  treatment  of  polyatomic  molecular  dynamics  an  entirely 
new  stochastic  theory  is  being  developed.  This  theory  is  designed  to 
handle  the  commonly  arising  cases  exhibiting  many  strongly  coupled  mole¬ 
cular  levels.  Both  intermolecular  and  intramolecular  situations  are 
being  treated.  In  addition,  a  vigorous  effort  is  under  way  to  quantum 
mechanically  handle  reactive  scattering  by  means  of  finite  element  theory 
formerly  developed  for  structural  mechanics  problems.  Other  specialized 
integral  equation  techniques  are  also  being  utilized  to  efficiently  solve 
the  dynamical  equations.  Finally,  a  general  sensitivity  theory  is  being 
pursued  to  provide  for  critical  parameter  analyses  in  macroscopic  kinetics 
as  well  as  in  microscopic  molecular  processes.  The  sensitivty  theory  is 
based  on  the  use  of  special  Green’s  function  methods,  and  the  approach 
shows  promise  for  application  to  a  number  of  engineering  and  chemical 
physics  problems. 

IV.  Research  Progress 

Outlined  below  are  the  specific  areas  of  progress  during  the  past 
year.  The  cited  references  are  from  the  list  of  publications  resulting 
from  this  contract. 

1.  Computational  kinetics  and  sensitivity  analysis  of  hydrogen- 
oxygen  combustion. ^ 

Kinetic  modelling  calculations  on  the  system  were  carried 

out  with  an  extensive  reaction  set  to  probe  the  vicinity  of  the 
three  explosion  limits.  Sensitivity  analysis  was  used  throughout 
this  investigation  tr  study  system  behavior,  in  particular,  to 
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elucidate  mechanistic  details.  The  concentration  and  sensitivity 
profiles  were  discussed  in  light  of  appropriate  experimental 
results  and  existing  theories  of  hydrogen  combustion.  The  results 
indicate  the  present  model  to  be  useful  over  a  wide  pressure- 
temperature  range.  The  reaction  set  is  also  used  to  probe  the 
sensitivities  for  an  experimental  study  designed  to  measure  rate 
constants  of  an  important  elementary  reaction,  H  +  0^  +  M  -*•  HO^  +  M, 
involved  in  this  system.  The  versatility  of  the  reaction  set  was 
also  demonstrated  by  a  study  of  a  related  chemical  reaction,  the 
decomposition  of  hydrogen  peroxide.  Finally,  the  prospects  for 
utilizing  the  methods  and  results  of  this  study  to  examine  other 
complex  kinetic  schemes  were  considered. 

2.  Chemical  kinetic  functional  sensitivity  analysis:  elementary 
2 

sensitivies. 

Sensitivity  analysis  was  considered  for  kinetics  problems 

defined  in  the  space-time  domain.  This  work  extended  earlier 

temporal  Green's  function  method  to  now handle  calculations  of 

6ui 

elementary  functional  sensitivies  —  where  u^  is  the  i-th  species 

concentration  and  cu  is  the  j-th  system  parameter.  The  system 

parameters  can  include  rate  constants,  diffusion  coefficients, 

initial  conditions,  boundary  conditions,  or  any  other  well-defined 

variables  in  the  kinetic  equations.  These  parameters  are  generally 

considered  to  be  functions  of  position  and/or  time  resulting  in 

6u  (x,t;x',t') 

functional  derivatives  — s —  which  measure  the  sensi- 

tivity  of  species  i  at  position  x  and  time  t  with  respect  to  a 
variation  in  the  j-th  parameter  at  position  x'  and  prior  time  t'. 
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Derivation  of  the  governing  equations  for  these  sensitivities 
and  the  Green's  function  was  presented.  The  physical  interpreta¬ 
tion  of  the  Green's  function  matrix  elements  was  elaborated.  A 
given  element,  G^ (x,t;x' ,t')  can  be  interpreted  as  a  generalized 
memory  function  which  has  important  physical  consequences  for  the 
design  of  kinetic  systems.  A  family  of  reduced  sensitivity  co¬ 
efficients  was  defined  in  terms  of  integrals  over  the  functional 
sensitivities  such  that  a  variety  of  laboratory  questions  can  be 
addressed.  Finally,  a  discussion  of  the  relation  of  this  work  to 
earlier  research  was  presented. 

3 l  The  selective  preparation  of  excited  vibrational  states  using 

3 

the  stimulated  resonance  Raman  effect. 

A  non-perturbative  treatment  of  stimulated  resonance  Raman 
(SRR)  scattering  was  developed  in  order  to  study  the  feasibility 
of  laboratory  experiments  for  the  selective  preparation  of  highly 
excited  vibrational  molecules.  The  theory  is  based  on  the  deri¬ 
vation  of  an  effective  Hamiltonian  for  the  interaction  of  the 
electric  field  with  molecular  electronic-vibrational-rotational 
states.  The  resulting  reduction  in  the  number  of  coupled  differen¬ 
tial  equations  describing  the  SRR  effect  allowed  for  a  convenient 
numerical  solution,  which  is  important  in  intense  laser  fields. 

In  addition,  when  the  spectral  transitions  are  power-broadened 
to  a  considerable  extent,  the  solutions  become  analytic.  An 
application  of  the  theory  to  the  I-  system  shows  the  possibility 
of  selectively  populating  high  vibrational  levels  (v*  >  10)  in  the 
ground  electronic  state.  A  similar  analysis  on  the  O2  system 
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shows  favorable  laboratory  conditions  for  pumping  states  v'  =  4~12. 
Most  importantly,  it  now  appears  possible  to  selectively  populate 
wide  choices  of  vibrational  levels  by  suitable  variations  in  the 
frequencies  and  proper  selection  of  the  intensities  and  durations 
of  the  applied  laser  fields.  Active  design  of  laboratory  appara¬ 
tus  is  underway  in  conjunction  with  this  study. 


4.  Further  developments  and  applications  of  sensitivity  analysis 

4 

to  collisional  energy  transfer. 

This  work  considered  the  sensitivity  of  integrated  collision 
cross  sections  to  the  variations  of  the  intermolecular  potential 
parameters.  The  study  was  concerned  with  the  scattering  of  an 
atom  with  a  linear  rigid  rotor  and  an  atom  with  a  breathing  sphere. 
Attention  was  focused  upon  first  order  sensitivity  coefficients 
(i.e.,  the  gradient  of  the  cross  sections  with  respect  to  potential 
parameters) ,  from  which  an  entire  family  of  derived  sensitivity 
coefficients  may  be  obtained.  The  original  elementary  set  of 

where  ou  is  the  i-th  system 


sensitivities  are  expressed  as 


30 , 
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parameter.  Within  the  family  of  derived  sensitivities  a  special 

3a 

class  of  coefficients  ■■■  -g—  was  shown  to  be  particularly  im- 

P 

portant  in  determining  the  extent  to  which  a  set  of  measurements 
is  able  to  define  the  parameters  of  an  assumed  potential.  In 

3o_mi 


addition,  the  derived  sensitivities 


3aA^' 


have  an  immediate 


bearing  ou  the  interdependence  for  different  possible  laboratory 
measurements.  Finally,  the  global  behavior  of  cross  sections  in 
parameter  space  was  examined  and  a  non-linear  interpolation 
formula  was  derived  utilizing  sensitivity  information.  The 
utility  of  this  latter  formula  is  illustrated  in  Figure  1. 


5.  Sensitivity  analysis  of  rotational  energy  transfer  processes 
to  the  intermolecular  potential. 

This  research  considered  the  sensitivity  of  rotational  energy 
transfer  processes  to  variations  of  parameters  within  an  assumed 
model  intermolecular  potential.  The  following  cross  sections  were 
examined:  integral  state-to-state,  pressure  broadening,  effective 
diffusion  and  viscosity,  and  final  state  summed  integral  cross 
sections.  In  order  to  simplify  the  calculation  of  cross  sections, 
attention  was  restricted  to  the  scattering  of  an  atom  and  a  linear 
rigid  rotor.  Furthermore,  the  collision  dynamics  are  approximated 
by  using  the  infinite  order  sudden  (IOS)  method.  Particular 
attention  was  focused  upon  the  sensitivities  of  different  cross 
sections  and  combinations  of  cross  sections  to  the  various  para¬ 
meters.  The  first  order  sensitivities  were  also  used  to  derive 
new  coefficients  which  describe  how  the  potential  parameters  cor¬ 
relate  given  a  limited  set  of  cross  section  data.  These  coeffi¬ 
cients  are  shown  to  be  especially  important  in  determining  the 
degree  to  which  a  set  of  measurements  is  able  to  define  various 
parameters  of  an  assumed  potential. 

6.  Sensitivity  analysis  of  differential  cross  sections  to  the 
intermolecular  potential.^ 

This  research  treated  the  sensitivity  of  both  final  state 
summed  and  state-to-state  differential  cross  sections  to  variations 


of  parameters  within  a  model  intermolecular  potential.  In  order 
to  simplify  the  calculations  of  the  cross  sections  and  first  order 
sensitivty  coefficients  (i.e.,  the  partial  derivative  of  the  cross 


sections  with  respect  to  potential  parameters  (attention  was 
restricted  to  the  scattering  of  an  atom  and  a  rigid  rotor.  Figure 
2  gives  an  example  of  a  particular  cross  section  and  its  associated 
sensitivity  derivatives.  Special  emphasis  was  given  to  an 
examination  of  the  sensitivity  of  angular  features  of  the  cross 
sections  to  potential  parameter  variations.  This  is  facilitated 
by  fitting  the  cross  sections  to  functional  forms  which  contain 
several  adjustable  parameters,  each  of  which  controls  a  particular 
feature.  The  first  order  sensitivity  coefficients  were  then  used 
to  derive  quantities  which  measure  the  sensitivity  of  a  feature 
parameter  with  respect  to  variations  of  a  potential  parameter.  The 
first  order  sensitivities  were  also  used  to  obtain  so-called 
derived  sensitivity  coefficients;  this  particularly  involved  those 
which  describe  how  the  potential  parameters  are  interrelated  for  a 
given  set  of  differential  cross  section  measurements.  The  behavior 
of  these  coefficients  was  examined  as  the  angular  range  of  the 
measurements  was  varied  and  as  account  was  taken  of  the  finite 
angular  resolving  power  of  molecular  beam  detectors.  The  results 
were  used  to  determine  the  degree  to  which  a  set  of  differential 
cross  section  measurements  is  able  to  define  various  parameters 
of  an  assumed  potential. 

7.  Action  angle  variables  for  the  quantum  three  particle  system. ^ 
This  work  used  the  newly  developed  quantum  action  angle  variable 
formalism  to  reduce  the  three-atom  interaction  system  to  a  minimum 
of  four  independent  coordinates.  A  transformation  was  also 
presented  which  converted  the  Hamiltonian  to  a  form  amenable  to 
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the  direct  numerical  integration  of  the  resultant  Schrodinger 
equation.  Computational  considerations  relevant  to  the  application 
of  finite  elements  were  also  discussed  and  illustrated  with  a 
numerical  calculations  of  the  quantum  energy  levels  of  a  rigid 
asymmetric  rotor.  The  ordinary  differential  equation  derived  for 
this  last  problem  permits  an  unambiguous  assignment  of  the  quantum 
expression  for  the  total  angular  momentum. 

8.  The  collision  of  two  linear  rotors:  a  scaling  theoretical 

g 

analysis  of  the  H^-H^  anc*  systems . 

An  energy  corrected  sudden  scaling  theory  was  applied  to  rota¬ 
tion-translation  and  rotation-rotation  energy  transfer  in  the 
and  HF-HF  systems.  The  scaling  predictions  are  in  very  good  agree¬ 
ment  with  the  exact  quantal  values.  Two  physically  important  results 
are  particularly  transparent  in  the  scaling  analysis:  (1)  The 
well-known  effectiveness  of  ortho-H^  in  rotational  relaxation  can 
be  traced  directly  to  the  higher  tensor  order  interactions  which 
are  also  operative  in  para-l^  for  j  >  0;  and  (2)  transitions  in 
which  rotational  quanta  are  conserved  differ  appreciably  from 
those  in  which  rotational  quanta  are  not  conserved.  Scaling 
predictions  using  different  values  for  the  average  collision  range 
£c,  were  compared  to  the  exact  quantal  cross  sections.  The  results 
indicated  that  essentially  the  same  value  of  will  be  determined 
by  inspection  of  the  dynamical  solutions  and  by  the  scaling  analysis. 
This  point  is  of  importance  for  future  applications  of  the  scaling 
theory  to  the  inversion  of  laboratory  data. 


9.  A  comparison  between  finite  element  methods  and  spectral 

9 

methods  as  applied  to  bound  state  problems. 

The  finite  element  and  spectral  methods  were  applied  to  two- 
dimensional  bound  state  problems.  A  comparison  of  the  spectral 
method,  which  requires  a  global  basis  set  expansion  of  the  wave- 
functions,  and  the  finite  element  method,  which  requires  no  such 
expansion,  was  performed.  A  procedure  was  given  for  formulating 
the  finite  element  approach  and  for  achieving  fast  and  accurate 
results.  The  convergence  of  the  finite  element  calculations  was 
considered  and  shown  to  be  well  behaved.  The  basic  conclusion  of 
this  work  is  that  the  finite  element  and  spectral  methods  are 
complementary.  In  the  regime  where  an  accurate  zero-th  order 
Hamiltonian  can  be  identified  the  spectral  method  excels.  On  the 
other  hand,  finite  element  theory  is  better  for  complex  potentials 
that  do  not  lend  themselves  to  conventional  eigenfunction  expansions. 
A  discussion  of  the  extension  of  the  finite  element  method  to 
higher  dimensions  was  also  included  in  this  work. 

10.  An  improved  phase  shift  approach  to  the  energy  correction 
of  the  infinite  order  sudden  approximation.^ 

A  new  method  was  developed  for  obtaining  energy  corrections  to 
the  infinite  order  sudden  (IOS)  approximation  by  incorporating  the 
effect  of  the  internal  molecular  Hamiltonian  into  the  IOS  wave  func¬ 
tion.  This  was  done  by  utilizing  the  JWKB  approximation  to  transform 
the  Schrodinger  equation  into  a  differential  equation  for  the  phase. 
It  was  found  that  the  internal  Hamiltonian  generates  an  effective 


potential  from  which  a  new  improved  phase  shift  is  obtained.  This 
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phase  shift  is  then  used  in  place  of  the  IOS  phase  shift  to  generate 
new  transition  probabilities.  As  an  illustration  the  resulting 
improved  phase  shift  (IPS)  method  was  applied  to  the  Secrest- 
Johnson  model  for  the  collinear  collision  of  an  atom  and  diamtom. 

In  the  vicinity  of  the  sudden  limit,  the  IPS  method  gives  results 
for  transition  probabilities,  significantly  better 

agreement  with  the  'exact*  close  coupling  calculations  than  the 
IOS  method,  particularly  for  large  An. 

11.  Stochastic  theory  of  intennolecular  energy  transfer  in 

in  the  presence  of  radiation. ^ 

A  theory  for  internal  energy  redistribution  in  polyatomic 
molecules  perturbed  by  strong  radiation  fields  was  constructed. 

Use  was  made  of  stochastic  theory,  which  assumes  a  random  phase 
approximation  is  valid  after  appropriate  time  intervals.  This 
approximation  permits  the  Schrodinger  equation  to  be  replaced 
with  a  finite-difference  master  equation  for  the  probabilities 
of  occuyping  the  various  quantum  levels.  The  semiclassical  theory 
of  the  radiation-molecule  interaction  was  employed  in  this  work. 

The  laser  linewidth  and  rotational  effects  are  incorporated  into 
the  formalism.  At  each  stochastic  step  the  energy  changes  in  the 
molecules  due  to  radiation  are  estimated.  Thus  energy  conservation 
is  explicitly  taken  into  account.  Model  calculations  for  SC^, 
whose  transition  dipole  moment  matrix  elements  and  anharmonic  force 
constants  have  been  determined,  indicate  a  complex  interplay  of 
anharmonic  and  radiative  coupling.  Power  density,  laser  linewidth 
and  detuning  from  resonance  are  observed  to  play  a  significant  role 


in  energy  redistribution 
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Figure  1.  An  illustration  of  collisional  parameter  scaling  for  two  breathing 
sphere  cross  sections.  The  parameter  6  controls  the  vibrational 
coupling  in  the  system.  The  solid  curve  is  obtained  by  laboriously 
repeating  the  calculations  at  new  6  values,  while  the  dashed  curve 
is  a  result  of  using  the  sensitivity  analysis  generated  scaling  theory. 
Excellent  results  are  obtained,  particularly  in  comparison  to  using 
perturbation  theory  indicated  by  the  dash-dot  line. 

Figure  2.  The  top  panel  shows  the  0+4  differential  cross  section  and  the 
panels  below  it  are  sensitivity  coefficients  with  respect  to  the 
overall  potential  strenth  e,  the  well  depth  anisotropy  a,  the 
range  parameter  Rm  and  finally  the  range  anisotropy  b.  The  minimum 
in  the  cross  section  near  80°  is  due  to  rotational  motion  inter¬ 
ference  and  the  striking  structure  in  the  sensitivities  at  this 
angle  can  be  readily  used  to  interpret  the  physical  origin  of  this 
structure. 
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